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Abstract 

j^ . The results are presented of a numerical simulation of the switching 

Ch ' magnetic junction by a spin-polarized current pulse under applied mag- 

netic field with the current density and field below the threshold values. 
r ^ • A possibility is shown of the switching with controllable delay relative to 

the current pulse. 

O' 

The switching magnetic junction by spin-polarized current [T] attracts 
continuous attention both of experimentalists and theorists. This is due 
^- ' to a possible practical application of the effect for high-density recording 

ON , information by current in magnetic media, as well as interesting problems 

related with the attempts of theoretical interpretation of the experimental 
results and prediction new effects. 

There are two important problems in experimental investigation and 
£ ■ application of the effect, namely, lowering the switching threshold and 

^^ ' shortening the current pulses causing the switching. Solution of the prob- 

^^ . lems will allow decreasing the Joule heating of the junction layers and 

avoiding destruction of the switching elements by the flowing current. 

Nowadays several ways have been proposed to lower the threshold cur- 
rent density, under which the initial configuration becomes unstable and 
switches to another one with different electric resistance. These are proper 
choice of the barrier layer material [2J, using diluted magnetic semicon- 
ductors as the pinned and free layers of the magnetic junctions 3 , using 
combination of the layer parameters corresponding to efficient spin injec- 
tion from the pinned layer to the free one and "locking up" the injection 
at the interface between the free layer and the nonmagnetic layer clos- 
ing the electric circuit [4], using a noncollinear initial configuration of the 
magnetic junction [S] HO [S] . 
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In Ref. [9] attention was paid to a possibility of substantial lowering 
of the current density needed to switch a magnetic junction under joint 
action of the current and magnetic field with the current density and field 
below the corresponding threshold values when the factors mentioned act 
singly. If the magnetic field strength is close enough to the threshold value 
(the anisotropy field), but slightly smaller than it, then a current with 
density below substantially than the threshold value without magnetic 
field can help to switching. Note, that the presence of the magnetic field 
does not break local character of the exchange switching the magnetic 
junction, because the magnetic field of such a value cannot do switching 
alone (without current). 

Together with the lowering current density of the stationary current, 
the presence of the magnetic field opens possibility of switching by short 
current pulses. Varying the pulse amplitude and duration allows switching 
with a time delay relative to the current pulse. 

Theoretical investigation of the switching by short current pulses is 
related with solving a set of nonlinear differential equations (the Landau- 
Lifshitz equations) with variable coefficients. Because of difficulty of this 
task, it is need to use numerical solution. 

Let us consider a magnetic junction with the free layer thickness small 
compared to the spin diffusion length, so that the macrospin approxima- 
tion [9] is applicable. The high spin injection conditions [4] are assumed 
to be valid, so that the switching magnetic junction is determined by cre- 
ating nonequilibrium spin polarization with injected spins |10l 1 1 1] . Un- 
der these approximations, the free layer magnetization in presence of the 
time-depended current j(t) perpendicular to the layers is described by the 
following equation [9]: 
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Here M = M/|M| is the unit vector along the free layer magnetization 
vector M, H is the external magnetic field, H^ is the demagnetization 
field, H a is the anisotropy field, n is the unit vector along the anisotropy 
field, Mi is the unit vector along the pinned layer magnetization, k is the 
Gilbert damping constant, 7 is the gyromagnetic ratio; 

eH a L 
JO = pp (2) 

has the meaning of the threshold current density for the switching mag- 
netic junction by a stationary current (an infinitely long pulse) without 
magnetic field; /j,b is the Bohr magneton, a is the sd exchange interac- 
tion constant, r is the electron spin relaxation time, L is the free layer 
thickness, Qi is the spin polarization of the pinned (injecting) layer con- 
duction. 

The last term in the left-hand side of Eq. |T|) describes effect of the 
spin-polarized current on the magnetic lattice under high spin injection. 



Let us consider a configuration with H = {0, 0, H}, n = {0, 0, 1}, 
H d = ~4tvM{M x , 0, 0}, Mi = {0, 0, 1}; x axis being directed along 
the current, yz plane being coincided with the layer plane. On the 
spherical coordinates (9, cf>) with the polar axis along z axis the vector 
equation JT} determining direction of the free layer magnetization vector 
M = {sin 9 cos <jf>, sinSsin^, cos 6*} takes the form of the following set of 
equations: 

|Uyf^{-^(0,<M) + B«>)}, (3) 

Zr = IT^ { ^' <M) + AcB(<m ' (4) 

where 

A(0, <f>, t) = h + h a cos + cos cos 2 <f> + h a ^- , (5) 

JO 

B(<f>) — cos <j> sin 4>, (6) 
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A numerical solution of Eqs. ((3j), was found using Simulink software 
in the MATLAB system [T5] . 

The following parameter values were given: n = 0.03, h a — 0.01. The 
initial relative orientation of the pinned and free layers was assumed to 
be antiparallel (9 — n). At high enough current density and/or magnetic 
field such a configuration becomes unstable, however, an initial deviation 
is needed which is promoted by thermal fluctuations in reality. The pres- 
ence of the thermal noise creating deviation of the free layer magnetization 
from the initial unstable stationary equilibrium state was imitated with 
giving a small initial deviation from such a state by angle of 0.01 radi- 
ans in the layer plane where the demagnetization field does not prevent 
fluctuation deviation, so that minimal fluctuation energy is needed. The 
time dependence of the spin polarized current density j(t) was given as a 
rectangular pulse of w duration. The chosen values of the magnetic field 
and the pulse amplitude, H = 0.9H a and j = 0.9?'o, respectively, being 
90% of the corresponding threshold values, avoid possibility of switching 
the junction by the field or by the current singly. Therefore, the switch- 
ing describing by the solution to be found is caused by joint action of the 
current and magnetic field. 

The simulation results are presented below as a time dependence of 
the free layer deviation angle from the initial antiparallel orientation. The 
dimensionless time is laid off as abscissa with to = (1 + k 2 )/{A-k^M) as 
the time unit; at M — 900 G one nanosecond corresponds to 200 scale 
divisions of the abscissa (to = 0.005 ns). The step shows the current 
turning on and turning off instants of time. 

There is a minimal pulse duration under which the switching can occur 
yet with the given values of the magnetic field and the pulse amplitude. In 
Fig. 1 the free layer deviation angle from the initial antiparallel orientation 
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Figure 1: Time dependence of the deviation angle of the free layer magnetization 
vector from the initial antiparallel orientation under the switching with joint 
action of the current pulse and the magnetic field, j/jo = 0.9, H/H a = 0.9, 
w = 67.15t . 
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Figure 2: Time dependence of the deviation angle of the free layer magnetization 
vector from the initial antiparallel orientation with the pulse duration below the 
threshold value. j/j = 0.9, H/H a = 0.9, w = 67.14£ . 



is shown as a function of time at w — 67.15to- It is seen that the switching 
takes place with substantial time delay relative to the pulse. At the given 
values of the magnetic field and the current pulse amplitude, that w value 
is threshold one: with decreasing the pulse duration only by O.Olto the 
switching does not occur (Fig. 2). Under increasing in the pulse duration 
the switching delay time decreases (Fig. 3). 

Lengthening the pulse makes it possible switching with lower values 
of the magnetic field and/or the current density. In the limiting case of 
the pulse duration much longer than the precession period, the switching 
condition takes the form corresponding to the stationary current [14] . 
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It is seen from Fig. 1 that the switching occurs as two stages, and 
the longer stage is realized after the current turns off. This fact may be 
explained by the following way. If the initial position of the magnetization 
vector is not exactly parallel to the external magnetic field, but is deviated 
from the antiparallel direction by some angle if), then the magnetic field 
H > H a cos ip is needed for the switching, the lower, the greater the 
deviation. The role of the current pulse is that it, acting together with 
the magnetic field, deviates the magnetization vector from the antiparallel 
direction by an angle large enough that the magnetic field could do the 
junction switching alone. 
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Figure 3: Decreasing (cf. Fig.[lJ time delay of the switching under pulse length- 
ening, j/jo = 0.9, H/H a = 0.9, w = 75t . 



The results show possibility of the switching magnetic junctions by 
joint action of the magnetic field and current pulse. Varying the pulse 
amplitude and duration allows obtaining, if necessary, a delayed switching 
with controllable delay time. 
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